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Preface 


This is the second quarterly report on Contract NAS5-21798, "Studies 
of Soundings and Imaging Measurements from Geostationary Satellites." 
During the past three months, good progress has been made on all tasks. 
The prototype Man-computer Interactive Data Access (McIDAS) equipment 
has operated well following repair of both the tape drives and the video 


disk recorder. 



Introduction 


This second quarterly report covers work performed from 1 November 1972 

through 31 January 1973. In summary, progress by task is as follows: 

Task A . Investigations of Meteorological Data Processing 
Techniques 

Significant progress has been made in analysing and 
differentiating the several forms of line start errors in 
ATS image data. An earth edge detection algorithm has been 
thoroughly tested and is producing good data for error 
analysis. Comparisons of earth edge errors and landmark 
displacements show very close correlation, which confirms line 
start jitter as a major error source in wind computation. Initial 
tests indicate line start errors probably must be removed 
before correlation lag computations are made, but the reasons 
for these findings are not fully understood at this time. Work 
during the next quarter will concentrate on development of error 
removal techniques. Work on bispectral cloud height determination 
has been temporarily suspended pending acquisition of data suitable 
for testing concepts reported earlier. 

Task B . Sun Glitter 

Work on this task has been limited to description of the geometry 
involved because the graduate student working on this task has 
shifted her efforts to a different task in response to guidance 
from her advisory committee. Work should recommence during the 


next quarter. 



2 


Task D . Cloud Growth Rate 

The early analyses in which correlations were sought relating 
cloud thickness to cloud brightness as seen in ATS images 
produced experimental results significantly different from 
theoretical results reported by other investigators. Work 
during this quarter has concentrated on exploring the hypothesis 
that the departure of observed data from theoretical predictions 
may be caused by differences in backscatter characteristics of 
ice particles and water droplets. Results from this investigation 
should be available during the next quarter. 

Task E . Comparative Studies in Satellite Stability 

Parts IV and V of the mathematical description of the stability 

model being developed by Mr. Das were attached to the last two 

monthly progress reports. Included in this report are the appendices 

to Mr. Das.' report for those who now hold a set of the five 

parts previously submitted. The entire report, as well as additional 

information, will appear in the final report for this contract. 

Good progress has been made in coding the stability model for 
the 1108 computer and should be nearly completed by the end of 
the next quarter. 

During this quarter the serious mechanical problems which had limited our 
use of the McIDAS equipment were corrected and the system is now working well. 
The problem of earth edge, or "line start," jitter in the ATS data is now 
the most serious problem affecting our ability to produce highly accurate 
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quantitative information from ATS data. We are confident this problem 
is tractable and our major effort is devoted to its solution. 

On 1 February we were requested by competent NASA authority to 
restructure our efforts under this contract to accommodate a new task 
and we have agreed to do so. The objective of the new task will be to 
explore and to develop in concept meteorological applications of a two- 
meter aperture optic earth observation system in earth synchronous orbit. 
This task relates very closely to Task F, which was deleted by NASA from 
our original proposal^ and we welcome the opportunity to pursue this study. 

Task Progress 

Task A . Investigations of Meteorological Data Processing Techniques 
1. The Problem: 

Precision display hard copy ATS images frequently show obvious line 
start timing errors which appear as irregular deviations of the earth edge 
from a smooth curve. The characteristics of this jitter are quite variable. 
Line to line jitter of four to six elements (quite prevalent), sudden shifts 
of 10 to 15 elements, and slow meandering oscillations of up to 20 elements 
peak to peak have been observed. In some cases, all three forms are present 
simultaneously and present serious difficulties to wind computations based 
on frame to frame cloud motions, since the observed East-West motion is 
strongly affected by line start deviations. 

To screen data for problems in line start timing and to understand 
and correct for its effects, an earth edge detection algorithm was developed 
and tested. The algorithm appears to be quite effective in detecting the 
earth edge consistently. However, successful utilization of the earth edge 
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information in correcting wind computation has not yet been achieved. Results 
to date and plans for future effort on this problem are discussed in the 
following. 

2. Earth Edge Detection Algorithm: 

The basic algorithm used in the earth edge detection program is a 
threshold detector with noise compensation. Starting at element N (a 
substantial distance to the left of the left earth edge) the average noise 
for the space view is obtained by averaging samples N to N + 50. The 
algorithm next looks along the scan line 14 elements ahead (to N + 64). 

If the samples from N + 64 to N + 68 all exceed the sum of the mean space 
noise and a pre-set threshold value, then N + 64 is returned as the position 
of the earth edge. If the condition is not met, N is increased by one 
element and the whole process is repeated. The 14 element buffer between 
space noise averaging and threshold testing is required to prevent averaging 
of signal with noise. The requirement that five successive samples exceed 
the offset threshold is needed to prevent triggering on isolated noise spikes. 
The use of a threshold offset by the mean space noise prevents line to line 
(or frame to frame) variations in noise from causing corresponding variations 
in the position of the detected earth edge. 

The first tests with the algorithm were made with data -from day 204, 1969, 
which showed considerable variations in noise amplitude and line start position 
as well. In addition to testing the algorithm's sensitivity to noise, a search 
was made for the optimum threshold. The optimum threshold should cause the 
algorithm to trigger at an altitude of 20 km to 30 km above the surface. 

At significantly lower altitudes variations in cloud cover near the limb would 
cause erroneous variations in the position of the detected earth edge. A very 
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low threshold would cause triggering at a high altitude well above the 
clouds, but would also be very susceptible to noise in the signal. 

Thresholds of 5, 10, 20, and 40 digital counts were tested (the ATS 
signal data ranges from 0 to 255 digital counts, the maximum signal roughly 
corresponding to the radiance of the brightest cloud) . The 5 threshold 
displayed a significant correlation with noise variations and a large 
scatter as well. The 10 and 20 thresholds produced edges that were not 
noise correlated. The 40 threshold was not noise correlated, but did show 
large variations in edge position which were inconsistent with both expected 
jitter behavior and with the results of 10 and 20 thresholds. The 40 
threshold thus appears to cause triggering on clouds. As a result of this 
preliminary testing, the 10 and 20 thresholds were chosen for standard usage. 

In normal operation, the detection program finds edges for both 10 and 
20 thresholds. The mean of these two is used for the final detected earth 
edge and the difference is used to check the quality of the detection. 

The difference is typically about three samples + 2 RMS. 

3. Comparison of Edge Detection and Landmark Measurements: 

Landmark measurements made for navigation of ATS-3 images for day 204, 1969 
are uniquely suited for testing the consistency of the earth edge detection 
algorithm. The landmark used was a small island off the coast of South America, 
which occupied approximately one ATS field of view and which could be isolated 
to one line and within one or two samples (most landmark measurements are based 
on patterns extending over many lines and elements) . Neglecting earth rotation 
within the ATS data frame (this is justifiable for the case in question because 
the spin axis inclination is small) , the landmark element position and the 
earth edge element position, for the same line on which the landmark falls, should 
differ by a constant. For a five frame sequence, this was found to be true 
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to within about 1.3 elements RMS. This sequence of comparisons revealed 
that the anomolous landmark motion previously discovered was indeed a 
result of line start timing errors. 

4. Frame to Frame Comparison of Earth Edges: 

Comparison of detected earth edge position as a function of line number 
with a smooth earth edge calculated by the SSEC ATS navigation program 
reveal moderate amplitude (4-8 elements) line to line jitter, sudden large 
shifts (10-15 elements) taking place within a few lines and lasting for many 
lines; and, perhaps most disturbing, large amplitude (up to 20 elements peak- 
to-peak) meandering oscillations with a varying period of the order of 
several hundred lines. The last effect is not readily observable in hard 
copy images, but is extremely important in the computation of frame to 
frame cloud motions. . 

Since the earth edge irregularities (and the inferred line start 
irregularities) are not the same in each frame, part of the observed cloud 
motion in the data frame is due to differences in line start timing errors. 

For many of the day 204, 1969 images, this motion is so large as to be 
readily observable on the SSEC McIDAS prototype TV display system when a 
sequence of images is viewed as a movie loop. The general characteristics 
of this motion were found to be in agreement with that predicted from the 
earth edge behavior, i.e., the latitude regions and sign of the erroneous 
motions appeared to be consistent with the relative motion of the detected 
earth edge. Certain regions of latitude which had the same line start 
timing error should appear as nulls in the anomolous motion. The position 
of these nulls also appeared to agree with predictions. All these approximate 
visual checks indicated that an accurate numerical check with analysis of 
digital data was warranted. 
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5. Correction of WINDCO Lags : 

The WINDCO (wind computation) program uses either cross-correlation 
or Euclidean norm methods to measure the relative line and element motion 
(or lag) of a cloud mass from one frame to another. Assuming constant 
cloud velocity and perfect line start timing, the element lag between the 
first two frames of a three frame sequence (equally spaced in time) should 
be equal (approximately) to the element lag between the second two frames. 

When line start errors are present, the two lag results will be in 
disagreement. If the earth edge detection is consistent and the correlations 
are high, the disagreement in the two lags should be predictable from relative 
differences in earth edge position from frame to frame. 

To test this hypothesis, approximately 20 cloud tracers were selected 
and lag computations processed. If the lag discrepancies agreed with the 
line start difference discrepancies, the feasibility of correcting lag 
computations for line start timing errors would be demonstrated. Results 
to date have, in fact, not demonstrated this. The results appear to indicate 
that the opposite is true, although a number of error sources have not been 
fully investigated at this time, among which are choice of poor cloud tracers, 
false correlations, and cloud size effects. It may turn out that the line 
start errors have to be corrected in the data before the lag computations 
are made instead of correcting them in the lags afterwards. 
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Task B . Sun Glitter 

The following develops the geometry of polar satellite sun glint 
observations. The object is to find the location of the sun glint with 
respect to the satellite reference frame. 

Consider an earth centered reference frame; z toward the north pole and 
x and y in the equatorial plane. 

Define; point S location of sun 

point A location of satellite 

point P location of reflection 

R radius of earth 

R + h distance from center of earth to A 

point 0 center of earth 

w angle of incidence = angle of reflection 

I angle between earths equatorial plane and the 

plane of the orbit 

It is assumed that the positions of the sun and the satellite are known at 
all times, i.e., S(0 (x), <f> (t) and A(0 (x), cf> (t ) ) are known. 

o u A A. 

It was necessary to find P = P [ @ s (x)<|>g (x) 0^(x) $^(x) ] . 

If there were no wind the ocean would be smooth and the vector P would 

be in the plane described by S and A. Therefore, a relation between <|> 

and 0 is available. Furthermore, the normal to the surface at point P is 
P 

c A-P 1/2 

n = -j— -j t- _ where In] = (2-2cos2co) 

11 l A - y l 
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»!# 


Third, 


n 


1R 


A-P n 

= COSO) = 


A-P 


' RT 


so there are three quations and three unknowns, 0 ,.<f> , co. With this infor- 

P P 

mation, it is clear that the vector P can be found such that 

p . rA-P(0 a ,* a , 9 s .* s ) 


Now consider the satellite reference frame such that 


x' along heading line 

y* _L heading line 

z* toward the center of the earth 


Let the rotation matrix which rotates the earth centered reference frame into 
the satellite reference frame be represented by 0(t). Let a be any vector 
in the earth system and a' be the vector in the satellite system. Then 
a = 0(t)[a-A] = 0(t)a - (R+h)z'« With this matrix P', the position of 
the sun glint in the satellite coordinate system, can be calculated at all 
times. 

For easiest analysis of the sea slope the camera should be pointed directly 
along P, For this orientation the relationship between the photo coordinates 
and the isolines of constant slope is the simplest. 


cos0 TT sin0 .+sin9„sin8 ,cos0 . 
HA H A A 


-cos8 t ,cos 9 ,+sin0„sin0 , sin0 , -sin0, cos© . , 
H A H A A HA' 


Q(t) = \ -sin0. sin8 TT +cos8 Tl sin8 . cos8 . 

A H H A A 


sin0 cos 0 +cos0 TT sin0 . sin0 4 -cos0 TT cos0 . 
H A H A A HA 


cos8 ,cos0 , 
A A 


cos0 , sin8 . 
A A 


sin0 , 


cosl 

where cos8 = cos 0 and tan 0 

HA H 


/sin 2 I - sin 2 0, 


cosl 


®H is the angle between the heading line and the plane parallel to 
the equatorial surface. 
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Task D . Cloud Growth Rate Study 

To parameterize convection, one of the basic quantities which must 
be measured is the volume, flux through a cloud. Previous work on 
evaluating convective heat transport from a geosynchronous satellite, such 
as done by Sikdar (1971) , determined the volume flux by measuring the 
horizontal divergence of the cirrus outflow of convective towers. However, 
the thickness of cirrus anvil was not known so an assumed value was used. 

To improve the estimates of volume flux, the thickness of the cirrus outflow 
must be determined for each cloud measured. 

A satellite sees a strong brightness gradient on convective clouds. 

One parameter which determines the brightness of a cloud is its thickness. 

The theoretical work of E. Hewson as presented in Kondratyev's "Actinometry" 
shows the dependence of the cloud's reflectivity on cloud brightness. However, 
the reflectivity reaches 80% for clouds several hundred meters thick and 
remains virtually constant for cloud thickness over 600 m. Hence, thick clouds 
should be uniformly bright, but satellite's do not see this happening. The 
difference could lie in the fact that almost all theoretical work on scattering 
has been done using the sphere as the geometry of the scattering particle. At 
present no measurements of the volume phase matrix for ice crystal clouds has 
been made available ("Major Problems in Atmospheric Radition" - Atmospheric 
Radiation Working Group, Bulletin of the AMS October, 1972). If it were, 
present methods could provide solutions to radiative transfer problems. 
Experimental work, such as by Huffman and Thursby (JAS, 1969) has shown that 
ice crystals have different scattering properties than spheres. A theoretical 
study- by Liou (JAS, 1972) .using long horizontal randomly oriented cylinders 
has shown that cylinders have considerable less back scatter than polydispersed 


spheres . 
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S-ince a satellite is looking at the top of a cloud, the scattering 
particles which the satellite sees are mainly ice crystals for a deep cloud. 

If ice crystals and rain drops have different scattering properties, then 
the brightness gradient could be caused by the different amounts of the two 
water forms in a cloud. If this were so, then the angular scattering dependence 
of the bright parts of the cloud should be different from that of the rest 
of the cloud. Analysis of this problem is presently being worked upon. 

Analysis routines have been written and a test set of convecting clouds in 
the Caribbean has been selected. The McIDAS system has been used to obtain access 
to the data. Analysis is presently being performed on the data, but no results 
have been obtained yet. 

Task E . Comparative Studies in Satellite Stability 

The mathematical development of the stability model has been completed 
and the model is now being coded for the 1108 computer. The appendices to 
the previously submitted mathematical analysis report are attached to this 
progress report. 

New Technology 

No items of new technology have been produced during this quarter. 

Program for Next Reporting Interval 

Task A . Investigations of Meteorological Data Processing Technique 

Continued work on the problem of ATS line start jitter is directed at 
the following tasks: 

1) Theoretical calculations of atmospheric limb radiance profiles, 
including investigations of the effects of sun angle, latitude, 
ozone, aerosols and atmospheric state; 
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2) Improvement of the edge detection algorithm to take account 
of the results of task (1) and to be more efficient in 
terms of cost and time; 

3) Application of the algorithm to study line start error 
characteristics for a wider selection of ATS data; 

4) Futher effort to correlate east-west residuals in wind 
computations with line start errors inferred from earth 
edge detection; 

5) Development of an algorithm to detect the position of the 
right hand earth edge; and 

6) Implementation of earth edge detection in routine data 
processing. 

Task B . Sun Glitter 

Work on this task depends upon availability of qualified graduate student 
assistance. This task will be curtailed to provide funds for the new Task H. 

Task D . Cloud Growth Rate 

Attempts will be made to differentiate and measure relative backscatter 
characteristics of liquid and ice cloud forms using McIDAS. 

Task E . Comparative Studies in Satellite Stability 

The next quarter will be devoted entirely to coding and debugging the 
stability model in the 1108 computer. 

i 

Task H . Large Optics Applications Study 

A broadly based study effort will be organized and pursued within SSEC, 
in a graduate seminar, and in conversations with persons in NASA and other 
government agencies. 

Conclusions and Recommendations 

None at this time. 



APPENDIX 1 


The Form of and £ B 

* 

The (110 x 1) generalized position vector q , i = 1 - 110, 

consists of the time-dependent coefficients introduced with the assumed 

displacement functions of the structural elements of the body B. The 

* 

coefficients for each element of £ are listed below. 


Element Element 


Element 


Element 


No . Name 

No. Nj 

1 y 17,l 

19 3 

2 y 17,2 

20 

3 y 17,3 • 

21 1 

4 y 18,l 

22 2 

5 y 18,2 

23 3 

6 y 18,3 

24 1 

7 y 19 , 1 

25 2 

8 y 19,2 

26 3 

9 y 19,3 

27 

10 y 20,l 

28 1 

11 y 20,2 

29 2 

12 y 20,3 

30 3 

13 q l,l 

31 1 

14 l q l,2 

32 2 

15 2 q l,2 

33 

J 3 

16 3 q l,2 

34 

17 l q l,3 

35 1 

18 201,3 

36 2 


A-1.2 


Element 

No. 

Element 

Name 

Element 

No. 

Element 

Name 

37 

3 q 4, 2 

63 

X 2 ,7 

38 

l q 4,3 

64 

X 

N3 

oo 

39 

2 q 4,3 

65 

X 2 ,9 

40 

3 q 4,3 

66 

x 2,10 

41 

y 13 

67 

X 2,ll 

42 

y l4 

68 

X 2,12 

43 

y 15 

69 

X 2,13 

44 

y 16 

70 

X 2,14 

45 

y l 

71 

X 2,15 

46 

*1.1 

72 

X 2 ,16 

47 

*1,2 

73 

. X 2,17 

48 

y 2 

74 

X 2,18 

49 

*2,1 

75 

a 2,13 

50 

^2,2 

76 

a 2,14 

51 

a 2,l 

77 

a 2 ,15 

52 

a 2,2 

78 

X 2 ,1 

53 

a 2,3 

79 

E l.l 

54 

a 2,4 

80 

a l,2 

55 

a 2,5 

81 

a l»3 

56 

a 2,6 

82 

a l,4 

57 

a 2,7 

83 

a l,5 

58 

a 2,8 

84 

a l,6 

59 

a 2,9 

85 

a l,7 

60 

a 2,10 

86 

a l,8 

61 

a 2,ll 

87 

a l,9 

62 

a 2,12 

88 

a l,10 




A-1.3 


89 

*1,11 

90 

a l,l2 

91 

X l,8 

92 

X l,9 

93 

X l,10 

94 

X l,ll 

95 

a l,13 

96 

a l,14 

97 

a l,15 

98 

X l,7 

99 

X l,l 


100 

b ll 

101 

b io 

102 

b 9 

103 

b 8 

104 

b 7 

105 

b 6 

106 

b 5 

107 

b 4 

108 

b 3 

109 

b 2 

110 

b l 


The constant matrix [G^, ] is defined by 
— B = ^ G B^% 

In this equation, C is a (3 x 1) vector and is the shift of the 

B 

center of mass of the satellite. Hence [G ] is a (3 x 110) matrix. 

B 

The nonzero elements of [G^] are listed below. All elements are to 
be divided by - 1 1 M |[^ . 


FIRST ROW 


Column 

No. 

Element 

Column 

No. 

Element 

1 

m 17 


4 

*18 

7 

“19 


10 

*20 

13 

b p T* 


21 

1 Pi ^ 

— b p 2 (e -1) - b p 2 -j 

22 

— D • H 2 
2 b p 2 * 


23 

bV 11 

27 

- bV* 


35 

1 P l^ 

b p 4 (e - 1)+ b p 4' 11 

36 

107 

2 b P 4 * 
2a B s P B**s 


37 

- bV* 


SECOND 

ROW 



2 

m 17 


5 

*18 

8 

*19 


11 

m 20 

14 

1 P l^ 

7“ b p i (e ' b p i’ 1 


15 

2 b p l 

16 

bV* 


20 

- b p 2* * 

28 

1 P l^ 

" b p 3 (e -1) + b p 3’ £ 


29 

- - D -l 2 
2 b P 3 

30 

105 

-b P 3‘" , 

2a B s p B*^s 


34 

b p 4* fc 


A-1.5 


THIRD ROW 


3 

m 17 

6 

m l8 - 

9 

*19 

12 

m 20 

17 

1 ^1^ 

r b p i (e ' 1} - bV* 

18 

— p - £ 2 

2 b P l 

19 

o • £ 
b p l 

24 

1 p l £ 

^b P 2 (e '^'bV 

25 

— 0 • 5. 2 

2b P 2 1 

26 

b P 2*^ 

31 

1 p l* 

b P 3 (e _1) “ b P 3' £ 

32 

— 0 • » 2 
2 b P 3 




P l £ 

33 

bV* 

38 

b p 4 (e _1) " b p 4 

39 

— D ' 

2 b P 4 * 

40 

b p A*^ 

41 

m 13 

42 

m l4 

43 

“15 

44 

m 16 

45 

"l 

48 

m 2 

63 

“l p P 2 4 

64 

2 

2 p M 2 B 

65 

Equal to element 64 

66 

Equal to element 63 

67 

Equal to 64 

68 

Equal to 63 

69 

Equal to 63 

70 

Equal to 64 

71 

Equal to 64 

72 

Equal to 64 

73 

Equal to 64 

74 

Equal to 64 

75 

*a 

76 

IT 2 

p p i 3 

10 B p 2 

77 

7 r 2 

- — p- 

78 

2 

p_ 

15 B p H 2 


B p 2 

91 

“ a 2 p p l 4 

92 

Equal to 91 

93 

Equal to 91 

94 

Equal to 91 

95 

o 2 

2u p, a_ a_ 

96 

7T 2 

- — — a„ p. 

p H l B 3 


10 B p 1 

97 

IT 2 

~ 15 p p l 3 B 

98 

71 2 

“ 8 p p l a B 

99 

tt p, a 2 




?1 B 




A. 1-6 


In the above, 

1 2ir 


a = J J L r 

x 0 0 


/ / [r 2 - 2a’r cos6 + (a') 2 ]log 


r ~ r 2 -2a'r cos6 + (a') 2 1 

l-2a'r cos0 + (a’) 2 r 2 


rd0dr 


r 2lr . 2 „ 2, „ / r 4 - 2ar cos8 + a 4- ^ 

a_ = J J (r -2ar cos© + a )log| 2~2 / rc * 6c ^ r 

0 0 \ 1 - 2ar cos0 + a r / 


a 3 = -7 / J 0 (x)xdx; W = ° 

A 0 0 


= 0.215922106463. 


The distance of the singular points from the center of the plate are 
denoted by a and a' . ! 
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APPENDIX 3 

Reduced Equations for Beams 

Let the (3 x 7) functional matrix [, K ] be defined by 

b i 


W ■ t b k i 1 (b q i ) 


(A3 .1) 


Then the elements of [^k^] are given by 


W - 


10 0 
P 1 S 

0 (e -1) s 
0 0 0 


0 

1 


0 

0 


p l s 


0 (e -1) 


0 

0 


0 

0 

1 


(A3. 2) 


We also have 




0 0 0 
0 10 
0 0 1 


(A3. 3) 


Then from (Eqs. (1.22) and (1.23), we get 

l 

(J> -J! l [ b K.]ds 


*b A il* 


A. . p l s 


( J ~ JL) J LK ] e 1 ds 

1 0 b x 
i 

j LK. ]sds 
i— 0 b 1 


b p i 


(A3. 4) 


We now introduce the (3 x 7) matrices [,K ] and [,K. _], of which the 

b iZ b id 


the nonzero elements are given by 

r b K i2 ] 2,2 = " ( b E i )( b I i,3 )( b T i )p l 6 *1. „ sin “i,2 

1 , z 


(A3. 5) 


3 P 1 S * 


f b K i2 3 3,5 ~ W^b I i,2 ) Vi )p l 6 *T. sin a i,3 (A3. 6) 

’ * i,3 ’ 


3 P 1 S 


[ b K i3 ] 2,2 = WVi.^l 6 (P 1 + K T i>2 Sin a l,2> 


(A3. 7) 



Using Eqs. (A3. 5) and (A3. 6), the (7 x 7) matrix [^A^l is. obtained 


as 


[ b A i2 J 


* } - J i> {2 h 0i t \ l W ds - / [ b K i2^ ds 


b p lU * /"[ b K.]e PlS ds - /\ K i2 3e Kl ds } (A3. 9) 


Pi s 


_ {2 b p i W B l [ b K i ]sds " / [ b K i 2 J sds} 

From Eqs. (A3. 7) and (A3. 8), the (7 x 7) matrix ^A^] is obtained 


as 


[ b A i3 J 


VA + "B > / W d “ - / 0 'b K i3 lds) 

(J + V / ds - [ [ b K i3 lePl ds} 


t b p l < % + V- / <b K i lsds * l I b K i3 lsds} 


(A3. 10) 


It is to be noted that, in the absence of external forces, 


GX 


IH 


-BC 


(A3. 11) 


Therefore from Eq. (1.103), we get 


®B = U V "l ' T 2 ) °B a B + ( Y T 7 * V°bS| 


- (t 3 £+ x 9 e; (A3. 12) 

Hence from Eqs. (1.22) and (1.23), the (7 x 1) vector ^b A i4^ 


obtained as 
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*bV ■ 


C b P i ) {(c B +u g) / R^ds - — ^— [y®] (t, 9 + T fl e)} 


n 


j v j— g— y 


^■ J l )( b'l )1( 'i + ^V?! e Pl ds - < ^,-‘ 1 ' 1 < T 3 i + T 9— ^ > 


. •* ** 


C P - ) { (u> + OJ )/ R.sds - 

b 1 E B 0 - 2||M||g 




( t 3 ®_ + T gi)> 


(A3. 13) 


In the above equation 


' R i,l + S ^ 


R. = J R 


R 


1.2 

1.3 


(A3. 14) 


J 


The forms of the (7 x 3) matrices [X^3, [ 3 and [X^3 will 
now be obtained. 


[^3 


^ (e 1 

p 2 -p 

W 4- l-il 


(A3. 15) 


[x 2 ] = 


+ — — (^ T - T ) } 

B lu.h y 7 O 


' pi Nib 

bl 1 P 1 1 B ilMl^ Y 7 

2 

b^i' 2 


V } 


(,p.) -V tu®]£2^ B + -^- 4 y - v> 


IlMllg 


(A3. 16) 
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APPENDIX 4 

Remarks on the Plate Functions 

Several forms of the plate function ^(z,w) have been shown on 

page 30, part I of this work. It will now be shown that the function 
* 

ip^(z,w) is the best choice. It is easily seen that the functions 

* * . * 

<i>^(z,w), i = 3-6, are included in and ip^ . Therefore, we 

* 

now examine ^p^(z,w). 

Let P be a lateral concentrated force applied by the plate on 

an attached system. Let v and v be the shear forces in the 

x y 

plate due to the deformation caused by the force P in the direction 
of the y- and x-axes, respectively. Then 

P = - [v dx + v x dy] (A4.1) 

c 

where C is a small circle of radius e with its center at the 

I 

point of application of the load P. 

From Eq. (A4.1), we get 


P = D V 2 ijj*]dx + [-^ V 2 (J> 2 Jdy} 


(A4.1a) 


where V 2 is the Laplacian in x and y, and D is the stiffness 
of the plate, or 


2 * 2 * 

d\ d\ 

p ■ iB % < £ 1 + & — - ldy - 


3z 9z 


9z3z 


or 



P = -96ttD( 1 - r^) (2 - 3rQ>e 2 rQ cos 2e Q 


(A4.2) 



Thus from Eq. (A4.2), it is seen that the force expression has a 
2 2 

zero at r^ = — , which is incompatible with the real situation. 

* * 

Now using ip^(z,w) instead of ip^(z,w) in Eq. (A4.1a), the ex- 
pression for P with 9 q = 0 is 

P = 24irDe 2 rQ(l - r 2 ) (A4.3; 1.79) 

Equation (A4.3) does not have a zero for 



A-5.1 


APPENDIX 5 

Reduced Equations for Plates 
a) Plate Nos. 1 and 4 

Substituting Eq. (1.34) in Eq. (1.33), we get 


+ W r)a i,13 + (r - 1)r2 S i,14 + Cr ' 1)r3 S i,15 


+ [^(A-jT)^ 7 + (r-l)r a ± g + (r-l)r £ ± g ]cos0 + [J^A.^)^ 


10 


+ (r-l)r a ± i;L + (r-l)r a ± 12 ]sin0 + [^(A^)^ 1 + (r-l)r 2 
+ (r-l)r 3 a ± 3 ]cos 26 + [J 2 (A 2 r)A‘ i 4 + (r-l)r 2 a^ 5 
+ (r-l)r 3 S 16 ]sin 20 + (r 2 log r)ji i>7 + * 1>8 Si i>8 + V 9 X i>9 

+ ^,10*1,10 + ^i,ll^i,ll } + ( pV "3 (1 + pV ( p K i,0 )sin p^i^ 1110 

3 B 

+ cos 6 ) (~£ - l){A 3 J 1 (A 0 r)a 1>13 + 9a. jl4 + (32r - 9)a i>5 

a 2 - - ■ 

+ [ “ ^(A^JsinQ + A 2 (i J 1 (A 1 r) - A 1 J Q (A ;L r) )cos " ^sin0 a^g 

A 2 

3 A 1 ?! 

+ [ (— - 8)sin0 + 8 cos0]a^ ^ + [ — — J 3 (A 3 r)cos0 + A^(— J 3 (A^r) 

3 3 

- A.J n (A r))sin0]a. in H — cos0 a. + [8 sin© - ( 8)cos0]a. 10 

10 1' i,10 r i,ll r i,12 

2 2 

2X 2 2X 2 3 

+ [— — J 2^2 r ^ sln 20 + (— ^— J 2 (X 2 r) - A^ (A 1 r))cos 20^ ^ 


+ 5(cos 20-2 sin 20)a. _ +[(24r-5)cos 20 - (24r-10)sin 20]a 

i , Z i»3 

2 2 
• 2X 2 2X 2 3 

+ [ — J 2 (A 2 r)cos 20 + (— — J 2 (X 2 r) - A^.^ (A-^r) ) sin 2033.^ ^ 
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+ 5(sin 28+2 cos 29)a. _ + [(24r-5)sin 28 + (24r-10)cos 28 ]a, - 

i>° 


11 


+ 7 X . 7 + I X. , 7 [r 2 - 2rr cos(6-8 ) + vU X [1 - 2rr cos(0-6 ) 

^ 3=8 3 3 3 3 3 

+ r 2 r 2 ] 2 (l-r 2 ){r 2 (l+2r 2 ) - rr . [cos (6-6 . ) + sin(8-8 )] 

J J J J 3 3 

- r\ . (2+r 2 ) [cos (8-8 . ) - sin(6-0.)] + r 4 r 2 [cos 2(8-8.)+ sin 2(8-0.)]}} 
3 3 3 3 3 3 # 3 


2 ,- 1 . 


“ p p i (a) B,l + “B,2 )x i “ 4 ( p D i ){A 0 J 0 ( V )a i,13 + r a i,l4 + (64 " r )a i,15 

a B 


+ X^JL (X r)cos8 a. r- cos9 a. 0 +— — cos6 a. Q + X J (X 1 r)sin0 a. 

XX X 1 • / / X • o 4 — 1 * / XXX X « 

r r 

15 


10 


X. Sine a + X sin0 a • + \ 4W )cos 26 a i,l " T COS 29 a i,2 

r r 


+ T cos 29 a i,3 + 4W> si ” 29 a i,4 - T sin 29 a i,S + T ain 26 a i,6 ! 
* <1 + “B,2 )R i + (i B,2 - “B.l“B,3 )a B r cos9 


( “B,1 + “B,2“B,3 )a B r Sln9 ' p Pit^3 ,(G A + U “B 


+ — ( 7 h " t 2 ,ig b1b + [ “b + “b“b + ( v T i ' t 2> 1g b% 


K 

1 


Mt 


(t 9 8* + t 3 8*)} 


(A5.1) 


b) Plate Nos. 2 and 3 

Equations for the plates Nos. 2 and 3 are similar to the above 
equation. The only difference is that these will involve twelve more 
functions of the form t^(z,w) introduced in Eq. (1.38). 



